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Adiabatic methods are potentially important for quantum information protocols because of their 
robustness against many sources of technical and fundamental noise. They are particularly useful 
for quantum transport, and in some cases elementary quantum gates. Here we explore the extension 
of a particular protocol, dark state adiabatic passage, where a spin state is transported across a 
branched network of initialised spins, comprising one ‘input’ spin, and multiple leaf spins. We find 
that maximal entanglement is generated in systems of spin-half particles, or where the system is 
limited to one excitation. 


I. INTRODUCTION 

Techniques of adiabatic passage are of interest from fundamental viewpoints, and also for their robustness against 
many sources of technical noise [T]. The canonical example is stimulated Raman adiabatic passage, STIRAP, [5], 
where an electron is transported between two meta-stable states, by applying two laser pulses that couple the meta¬ 
stable states to an excited state, or occasionally a continuum lai. By applying the laser pulses adiabatically in 
the so-called counter-intuitive ordering, where the unpopulated transition is coupled before the laser that addresses 
the particle to be transferred, the transported particle can move between the two states without ever populating the 
excited state. In this way the protocol shows robustness against spontaneous emission, and because the population at 
any time only goes like the ratio of the Rabi frequencies of the two fields, STIRAP is also robust against fluctuations 
in the total laser intensity when the fields come from the same source. 

A conceptual breakthrough occurred in 2002 when Brandes and Vorrath reported a method to use STIRAP to 
transfer electrons between wells in a double quantum dot. This was significant as it was the first time that the use of 
engineered, as opposed to naturally occurring, systems was considered. Following spatial STIRAP, schemes have been 
proposed where only the spatial tunnelling interaction is varied to effect the counter-intuitive pulse sequence, and these 
are sometimes termed coherent tunnelling adiabatic passage, CTAP. Eckert et al. considered atomic transport through 
optical potentials [6], Siewert and Brandes proposed transport in superconducting networks [7], and Greentree et al. 
proposed electronic transport in quantum-dot and phosphorus in silicon systems [8j. Later proposals included Bose- 
Einstein condensates [illD], optical waveguides mm, Bose-Hubbard systems m, sonic systems El, polarisation 
m, and spin chains [MU]. Spatial adiabatic passage in spin chains is usually termed dark state adiabatic passage, 
DSAP, and is the subject of this work. 

The combination of designed Hilbert spaces with adiabatic processes leads to interesting new protocols and oppor¬ 
tunities. In particular, adiabatic passage enables quantum gates via the method outlined by Unanyan, Shore and 
Bergman [TM2^ : interferometers [3S], interaction-free measurements [21], and robust splitters [25k27] . 

Here we consider the combination of DSAP with branched spin networks to explore how entanglement is generated 
between the leaf nodes, i.e. end points of the network. We also note that adiabatic passage on lattice networks has 
been considered [2H]. Our aim can be understood with respect to Fig. [^(a), which shows a simple network where all 
of the spins are initialised in their lowest spin state, except for the ‘initial’ spin, which is set to some particular state. 
We perform the counter-intuitive pulse sequence (described explicitly below) with an XX -|- YY -|- ZZ Hamiltonian, 
and explore how the initial spin state is distributed between the leaf nodes. In its simplest form, this case has been 
considered as a means to create a superposition of a single particle amongst the leaf nodes |21| 1221 ISO] , and this case is 
isomorphic to the spin-half case discussed below. We extend this work by exploring networks of higher-spin systems 
and multiple leaf nodes. We observe more complicated behaviour, although such behaviour does not appear to lead 
to entangled states superior to the spin-half case. 

This paper is organised as follows: We first describe the Hamiltonian for DSAP with high-spin systems. We then 
show results for spin-half, spin-one and spin-three half networks. In all cases we show the adiabatic evolution and 
quantify the final entanglement generated. Finally we compare the protocols and offer perspectives. 


II. HAMILTONIAN 

In its most general form, we consider the adiabatic transport of spin information from an initial spin, L, to some 
entangled state of final (or leaf) spins, Ri to Rn, via a middle spin, denoted M, which takes the analogous role 
of the unpopulated excited state in STIRAP. We assume a uniform, time invariant magnetic field B aligned with 
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FIG. 1: (a) Schematic showing a network for four spin-one particles, with starting node, L, central node M and two 
leaf nodes i?i and i? 2 - Spin-spin couplings are indicated by the solid lines connecting the spins, (b) Dark state 
adiabatic passage is achieved by utilising the counter-intuitive pulse sequence. To transport spin information from L 
to i?i and i ?2 requires 11 ^( 0 ) and ^ The exact shape of the pulse is of little 

importance provided that adiabaticity is preserved and here we have chosen squared sinusoidal coupling for 

simplicity. 


the Z direction which energetically separates manifolds of different numbers of excitations. We assume time-varying 
inter-site couplings between the spins and the middle spin only, where the L — M coupling strength is (t) and the 
M — Rj coupling strength is . The Hamiltonian for our network can be expressed as 


h = Y.bjI 




H.c. 


a(jfT 




Z jZ 
M^Ri 


( 1 ) 


where the index i ranges across all spins {L, M, i?i...i? 7 v}i and the index j ranges across all leaf spins Ri to 
Jz, and J~ are the Z-projection of the spin operator, spin raising and spin lowering operators respectively. The 
Hamiltonian preserves the number of excitations meaning that subspaces of different numbers of excitation can be 
treated separately in the adiabatic passage, a is the coefficient of the always-on ZZ coupling and is to prevent 
occupation of spin M, and operates in a similar fashion to the role of central state detuning in Ref. m- The value 
of a was set to 1 x for all systems. 

DSAP of the spin state from the initial to the final site(s), utilises the counter-intuitive pulse ordering, i.e. the VIr. 
coupling should be high before U,r. For our numerical results, for simplicity we have chosen: 


(^) — f^max sin 



^Rj (^) — f^max cos 



( 2 ) 


where fmax is the total time for the transport, and 0 < t < fmax is the time. This is shown in Fig. [^b). In keeping 
with other adiabatic passage protocols, the exact form of the D is relatively unimportant providing that the counter 
intuitive ordering is maintained, i.e. <C Q,R^(t) as t —)■ 0 , ^ (t) as t ^ tmax, and the variation in the 

ft satisfy the adiabaticity criterion, so that for every pair of eigenstates IV'), Itj)) with respective eigenenergies if|^) 
and F;| 0 ) 




( 3 ) 
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Although our calculations are performed explicitly as a function of time and without assuming adiabaticity, all of the 
results we show below are in the adiabatic limit. 

While the Hamiltonian and pulse ordering can allow for adiabatic passage of the state in site L to the recipient 
spins Rj, it is important to stress that DSAP is also dependent on the choice of the initial state of all spins other than 
L. For simplicity, we will assume a state with all spins M and Rj prepared in their lowest spin projection with respect 
to the z axis. Another enabling state has all spins in the highest spin projections and yields completely symmetric 
results. It is possible that there are other states that will permit DSAP, as was the case for direct transport along a 
one-dimensional chain HH], although we will not explore such possibilities here. 

We quantify the DSAP by two methods. First is to identify the final states for various network configurations. 
Secondly, we examine the entanglement that is generated between the Rj nodes. Because we are dealing with multi¬ 
partite systems of varying dimensionality, there is no unique way to quantify entanglement. For simplicity, we choose 
to use the entanglement of formation EH [31], which is an entanglement monotone for pure states where 0 corresponds 
to no entanglement and 1 to maximal entanglement for that system. The entanglement of formation for a state \ip) 
is defined 


= -Tr(pi^, logaPflJ = -Tr(pi^, log^PR^) (4) 

where pR^ is the partial trace of |'0)('!/'| over subsystem i ?2 and pR^ is the partial trace of |'0)('0| over subsystem i?i. 
In the configurations with more than two recipients we quantify the entanglement of formation as being between one 
arbitrary leaf spin and all other recipient leaf spins. 


III. SPIN-HALF NETWORKS 

The case of spin-half networks with one excitation is formally equivalent to the case of multi-recipient adiabatic 
passage [li|3D] and so we consider it here purely for the purposes of review and to compare with our other results. 

We denote the spin projections of the spins with respect to their eigenvalue. We define |1) as the state with 
spin projection parallel to the external magnetic field, i.e. with spin projection 7i/2; and |1) as the state antiparallel 
to the external field, with spin projection —Ti/2. The middle and leaf (right) spins are assumed initialised in |1), and 
we study the adiabatic passage for different initial (left) spin states. The evolution is indicated in Fig. The general 
form of the dark state, i.e. the eigenstate which defines the adiabatic passage, for n leaf nodes for the case considered 
above is 


\Do) 


nHfl|llll...l) - VLl (I1111...1) -H I1111...1) -H ■ ■ ■ -H |lll...l)) 

-I- 


(5) 


Equation ([^ shows that for the two leaf configuration, an initial excitation at L (Hi = 0) is adiabatically transported 
to an equally weighted, maximally entangled state of the Rj {^Ir = 0). As the number of leaves is increased, we see 
the formation of W-like states, which have a lower entanglement (as quantified by the entanglement of formation) 
with increasing n. 


IV. SPIN-ONE NETWORKS 

Spin one networks are more complicated than spin-half, and are no longer isomorphic to the simple case of direct 
particle transport, although in certain cases a spin-one chain can be treated as a form of particle transport |33[ I34j . 
Spin-one DSAP faithfully transmits the state of an arbitrary qutrit from one end of the chain to the single recipient 
spin at the other end of the chain [18] . However, as we will show, the branched geometry maps the doubly-excited 
state partially into product states of two sites with one excitation each, rather than solely an entangled state of both 
excitations in one or the other site, thereby failing to preserve the integrity of the qutrit. 

As before, we define our basis states with respect to the Z projections of the individual spins. The states are 
|1), |0), |I) corresponding to the spin projections h, 0 and —h. As before, we initialise spins M and Ri to |T) (antiparallel 
to the Z field). In the spin-one system, transport from L to Rj can be achieved for either one excitation (|0)) or two 
excitations (|1)). 

When limiting the system to one excitation, i.e. choosing the initial state of spin L to be |0), the transport achieved 
in all configurations is the equivalent of the spin-half case and the entanglement of formation for both spin-one and 
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FIG. 2: Population evolution for spin-half networks with (a) two leaf nodes, (b) three leaf nodes, and (c) four leaf 
nodes. The excitation begins at L i.e. in the state |ll...l) and moves to an equally weighted superposition of all of 
the leaf nodes, which is a W state across the spins Ri. All initial and final states as well as the quantified 
entanglement related to each protocol is summarised in Table 


spin-half is the same. The evolution for three configurations, with two, three and four leaf nodes, is shown in Figs.|^a), 
(b) and (c). 

In the two-excitation subspace for two leaf nodes, the evolution becomes more complicated due to the number of 
ways the excitations can be shared between the leaf nodes. Figure [^d) shows the transport with starting state |liil) 
with two leaf nodes. The counter intuitive pulse sequence ensures that all of the excitation is transported from L to 
the Rj without final occupation of M. However the hnal state is complicated by the fact that there are more ways for 
the excitation to be shared between the Rj than in the one excitation case. The pertinent eigenstate for the DSAP is 


I_ 2 H||iin) + (loiio) + |oM)) + ni ( 2 |iioo) + iiiii) + |ini)) 

° -1- -h 6^2 

This state shows features in common with more conventional adiabatic passage evolution. As 12^(0) ^ 12^(0), we 
retrieve the starting configuration, i.e. \Do) —>■ llllll). In the adiabatic limit, the state of M is always 1, and the 
final state has excitation only in the leaf nodes. The final state is, however, more complicated in the two-excitation 
subspace. We observe the formation of a superposition of a W-like state: |Illl) + and the state |1100), 

with most of the population in |1100). We find that |I2o) has less entanglement than a W state (as measured by 
entanglement of formation), and less entanglement than the state formed from the one-excitation state. 

The three and four leaf node configurations with two-excitations are more complicated, and we do not show the null 
states, as these are complicated and do not provide additional insight into the evolution. The final state conhguration 
for the three leaf node configuration gives rise to a superposition of two W-like states: 

|Z2o(tmax)) = ^ (inlOO) -f IITOIO) + IITOOI)) + ^ (lilin) + IlillT) + lITill)). (7) 

This entangled state includes all of the possible ways that the two excitations can be shared between the leaf nodes, 
but the system is more strongly weighted towards configurations where the excitations are shared most evenly. We 












5 


One Excitation 



(d) 


(e) 


a 


(f) 


;3 

a 

o 

Ch 


Two Excitations 



^/^max 


FIG. 3: Spin-one networks for DSAP. For the case of one excitation and (a) two leaf nodes, (b) three leaf nodes, and 

(c) four leaf nodes; we observe results equivalent to those seen for the spin-half networks. For two excitations and 

(d) two leaf nodes, (e) three leaf nodes, and (f) four leaf nodes, we observe evolution to a more complicated final 
state, via a transient, W-like state, as shown. The sharing of the population between the states is determined by the 
egalitarian principle as discussed in the text. Final states for two excitations and four nodes are not given here due 

to the large number of states involved, but is quantified, along with entanglement of formation, in Table 


term this sharing the egalitarian principle, i.e. we expect the adiabatic passage to always favour the configurations 
where the excitations are most equally shared between the leaf nodes. 

Similarly, the four leaf node configuration with two excitations also obeys the egalitarian principle, with a contri¬ 
bution from a W-like state, and states where the excitations are shared on two of the leaf nodes with an amplitude 
twice that of the W-like state. 


lA'oWmax) =^(1110011) -f IllOlOl) -f IllOllO) -f IlllOOl) -f IlllOlO) -b IllllOO)) 
V7 


-b (inilll) -b IlnilT) -b lITiliT) -b inilli)). 

2v 7 


( 8 ) 


As before, the system is weighted towards configurations where the excitations are shared as evenly as possible, and 
this results in a state with less entanglement than configurations with fewer leaf nodes. 


V. SPIN-THREE-HALF NETWORKS 

We now consider spin-three-half networks. We again define the eigenstates with respect to the Jz projection, 
labelling the states |3), |1), |1), |3) corresponding to the eigenvalues -b3?i/2, +h/2, —?i/2, —2,hl2 respectively. In the 
spin-three-half system, transport can be achieved by using one, two or three excitations. When limiting the system to 
one excitation, the transport achieved in all configurations is the equivalent of the spin-half case and the entanglement 
of formation was found to be equal in both cases. When limiting the system to two excitations. Figs. (a), (b) and 
(c), the transport produces final states that are equivalent to those obtained in the spin-one case. However the final 
states occur with different populations here than they do in the spin-one system. This is due to the additional -x/S 
factor from the and J~ spin operators, which also modifies the degree of entanglement produced. Such results are 
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FIG. 4: Evolution of systems of spin-three-half particles with two and three excitations. The two excitation results 
for (a) two leaf nodes, (b) three leaf nodes and (c) four leaf nodes, are qualitatively similar to the analogous spin-one 
systems, although quantitively different due to the different couplings. Three excitations for (d) two leaf nodes, (e) 
three leaf nodes, and (f) four leaf nodes, exhibit excitation sharing following the egalitarian principle as before, 
albeit with increased complexity due to the large number of states involved in the transport. The states typically 
not shown here for reasons of space, but all of the final states can be found in Table 


expected to quantitatively affect analogous transport through higher spin networks (i.e. networks with spin >3/2) 
without modifying the qualitative form of the adiabatic passage or the final entangled states. 

With three excitations (initial state |333 • • • 3)) the evolution is more complicated still, as shown in Figs. Kd), (e) 
and (f). However we still find that the egalitarian principle holds, with all possible final states represented in the 
final superposition, but the states with the most equal sharing of the excitation favoured. This sharing of excitations 
reduces the overall entanglement in systems with increasing number of leaf nodes. 


VI. DISCUSSION AND CONCLUSIONS 

Dark state adiabatic passage can be used to generate entanglement in small networks, and we have shown the 
properties of the evolution and entanglement in networks of spin-half, spin-one and spin-three-half systems. Our 
results for the final state and entanglement of formation for these networks are summarised in the Table |lj Networks 
of spin-half systems show the creation of W-like states across the leaf nodes, and can be understood analogously 
to particle splitting in spatial adiabatic passage networks [HJ |3S] . However the higher spin systems exhibit what 
we term an egalitarian principle, where the excitations are distributed across all possible leaf nodes, but with the 
states with the most equal distribution of the excitations more strongly represented in the final entangled state. 
This overrepresentation of states with equal sharing reduces the entanglement of formation observed in the high spin 
systems, with the consequence that systems with fewer excitations and fewer leaf nodes exhibit greater entanglement. 
This result is likely to be significant in the development of quantum protocols for high-spin, qudit systems. 
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TABLE I: Summary of final states following mnltiple recipient dark state adiabatic passage, and entanglement of 

formation. 


Figure Spin Nodes tfo 


Entanglement 

2a 

1/2 

2 

IlIII) 

;^(|IIIl) + |IIlI)) 

1 

2b 

1/2 

3 

IlIIII) 

^(IlIIII)+ |IIIlI) + |IIlII)) 

0.9183 

2c 

1/2 

4 

lllllll) 

1 (lllllll) + lllllll) + lllllll) + lllllll)) 

0.8113 

3a 

1 

2 

lolll) 

^(|II0I) + |III0)) 

1 

3b 

1 

3 

lolill) 

;^(|II0II) + |III0I) + |IIII0)) 

0.9183 

3c 

1 

4 

lOlllll) 

^ (lllolll) + IIIIOII) + IIIIIOl) + IIIIIIO)) 

0.8113 

3d 

1 

2 

IlIII) 

yi(iiioo))+yi(iiiii)+iiiii)) 

0.8072 

3e 

1 

3 

IlIIII) 

(IIIIOO) + IIIOlO) + IIIOOl)) 

+ yj(|IIlII) + |IIIlI) + |IIIIl)) 

0.7264 

3f 

1 

4 

lllllll) 

^(lliooil) + IllOlOl) + IllOllO) 
+|III00I) + IIIIOlO) + IIIIIOO)) 

+ y4 (lllllll) + lllllll) + lllllll) + lllllll)) 

0.5152 

- 

3/2 

2 

11333) 

^ (13313) + I333I)) 

1 

- 

3/2 

3 

113333) 

^ (133133) + 133313) + |33331)) 

0.9183 

- 

3/2 

4 

1133333) 

1 (1331333) + 1333133) + |333313) + |333331)) 

0.8113 

4a 

3/2 

2 

11333) 

y|(|33II)) + yi (|33l3) + 13331)) 

0.9021 

4b 

3/2 

3 

113333) 

1 (133113) + 133131) + 133311)) 

+^ (|33l33) + |333l3) + |3333l)) 

0.7080 

4c 

3/2 

4 

1133333) 

y^A(|33ii33) + 1331313) + |333311) 

+1331331) + 1333113) + |333131)) 

+ ^(1331333) + 1333133) + |333313) + |333331)) 

0.4997 

4d 

3/2 

2 

13333) 

^ (|33ll) + |33Il)) + ^ (13333) + |3333)) 

0.9763 

4e 

3/2 

3 

133333) 

^|33IlI) + ^(|333Il) + |33lI3) 

+|33l3I) + |33Il3) + |33I3l) + |333ll)) 

+ ^(133333) + 133333) + |33333)) 

0.6297 

4f 

3/2 

4 

1333333) 

^(1331113) + 1331131) + 

I3313II) + 1333111)) + ^(1331133) 

+|33l3I3) + |33l33I) + |33Il33) 

+ 331313) + 331331) + 333113) 

0.3889 


+1333131) + 1333113) + |333131) 
+I33331I) + 1333311)) + ^(|333333) 


+1333333) + 1333333) + |333333)) 
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